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ABSTRACT 

Determining the type of matter that is inside a neutron star (NS) has been a long-standing goal of astrophysics. 
Despite this, most of the NS equations of state (EOS) that predict maximum masses in the range 1 .4-2.8 Mq are 
still viable. Most of the precise NS mass measurements that have been made to date show values close to 
1.4Mo, but a reliable measurement of an over-massive NS would constrain the EOS possibilities. Here, we 
investigate how optical astrometry at the microarcsecond level can be used to map out the orbits of High-Mass 
X-ray Binaries (HMXBs), leading to tight constraints on NS masses. While previous studies by Unwin and 
co-workers and Tomsick and co-workers discuss the fact that the future Space Interferometry Mission should be 
capable of making such measurements, the current work describes detailed simulations for 6 HMXB systems, 
including predicted constraints on all orbital parameters. We find that the direct NS masses can be measured 
to an accuracy of '^2.5% (l-cr) in the best case (X Per), to ^6.5% for Vela X-1, and to ^10% for two other 
HMXBs. 

Subject headings: accretion, accretion disks — equation of state — astrometry — instrumentation: interfer- 
ometers — stars: neutron — stars: individual(X Per, Vela X-1, V725 Tau, GX 301-2, 
SAX J0635.2H-0533, V830 Cen) 
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1. EVITRODUCTION 

Neutron stars are found in a large variety of astrophysical 
settings, such as in accreting binaries, degenerate binaries, 
supernova remnants, and as isolated objects. They are one 
possible endpoint to the evolution of massive stars as well 
as being the locations of the highest magnetic field strengths 
and highest densities in the Universe. They are often found 
through their radio. X-ray, or gamma-ray pulsations, which 
can be on time scales as short as mill iseconds, and they can be 
accre t ion- or rotation-pow ered ( Tavlor. Ma nchester & Lvna 
Il993t iBildsten et all Il997t IXbdo et al. 200i). Despite at- 
tempts to determine their fundamental properties using nu- 
merous techniques, the form of matter that exists within a 
neutron star is still a mystery. 

The work that has been done to try to determine the neutron 
star composition has strong contributions from both theory 
and observation. For an assumed neutron star (NS) composi- 
tion, the theoretical pressure-density relationship (i.e., equa- 
tion of state, EOS) directly predict s a spe cific NS mass-radius 
relationship (Lattim er & PrakashI 1200 It) , so that measure- 
ments of mass and/or radius provide direct constraints on the 
NS EOS. It is currently even unknown whether NSs might be 
made of quark matter, which gives a radi cally different mass- 
radius relationship from normal matter jLattimer & PrakashI 
|2004). Most notably, each theoretical EOS has a maximum 
NS mass that can be supported (Lattimer & Prakash 2004), so 
that one measurement of an over-massive NS (~2M0 rather 
than the more canonical 1.4Mq) would severely constrain the 
possibilities for the NS EOS. 

Although accurate NS mass measurements have been made 
for NSs in binary systems, the very accurate measurements 
have found values close to 1.4M0 (iThorsett & Chakrabartyl 
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Il999t iLattimer & PrakashI l2007h . However, more recent ob- 
servations have shown evidence for over-massive NSs with 
best estimates in the L8-2.8 M0 range. Th ese have corne both 
from High-Mass X -ray Binaries (HMXBs. iBarziv et alj|2001l: 
IClark et aP l2002|) and observations of binary radio pulsars 
(iFreire et alJl2008l) . In addition, there is evidence from stud- 
ies of "type I" X-ray bursts, wh ich are caused by t hermonu- 
clear flashes on the NS surface dLewin. van Paradijs & TaarnI 
1 19931). that the maximum NS mas s is in the 1 .9-2.3 M© range 
(ISteiner. Lattimer & Brownll2010l) . Improving the mass con- 
straints and enlarging the sample of NSs with accurate mass 
measurements could finally lead to definitively ruling out 
EOSs. 

Currently, the HMXBs with the best NS mass estimates use 
X-ray pulsations to measure the projected size of the NS or- 
bit and optical spectroscopy to measure the radial velocities 
of the companion star and, thus, to constrain the projected 
size of the companion's orbit (Charles & Coe 2006). How- 
ever, even in cases where these measurements are possible, 
the binary inclination is still a major source of uncertainty 
in the measurement of the NS mass. Thus, in this work, we 
explore how NS mass measurements can be improved via op- 
tical astrometry. In many cases, HMXB orbits are as large 
as 20 to several hundred /j,as, and the companion's orbit, in- 
cluding the binary inclination, can be m apped directly with 
astrometry on the microarcsec ond level (lUnwin et al J 120081 : 
'Toms ick. Shaklan & Panll2009h . Here, we look at the spe- 
cific case of the future Space Interferometry Mission and per- 
form detailed simulations to determine how well this mission 
would be able to constrain NS masses. 

2. SELECTING SOURCES FOR NEUTRON STAR MASS 
MEASUREMENTS 

The best targets for obtaining neutron star mass measure- 
ments by mapping out binary orbits with optical astrometry 
should meet several criteria. They should be optically bright, 
and the angular sizes of their orbits should be large. In ad- 
dition to the orbit being large, the cleanest orbital measure- 
ments will be obtained in cases where all or nearly all of the 
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Table 1 

Estimates of System Parameters and Astrometric Signatures 







d' 


Pj' 




Mns" 


P 


PcLimp 


Source Name 


V 


(kpc) 


(days) 


(Mq) 


(Mq) 


(pas) 


(pas) 


4U 0352+30/X Per 


6.6 


0.7-1.3 (1.0) 


250.3 


15.5^ 


1.4 


1993 


165.1 


3A 0535+262/V725 Tau 


9.6 


2.0"° 


111 


8-22" (15) 


1.4 


574 


49.0 


XTE J 1946+274 


16.9 


9.5 


169.2 


10-16"^'" (13) 


1.4 


153 


14.9 


Vela X-l/GP Vel 


6.9 


1.9 


8.96 


23.8' 


1.86-2.27'" (2.0) 


131 


10.2 


GX 301-2/BP Ciu 


10.8 


3.0^.0*'' (3.5) 


41.59 


43'° 


1.85' 


238 


9.8 


2S 1417-624 


17.2 


1.4-1 1. 1 (6.0) 


42.12 


12" 


1.4 


93.7 


9.8 


EXO 2030+375/V2246 Cyg 


19.7 


7.1 


46.02 


17.5'" 


1.4 


94.2 


7.0 


SAX J0635.2+0533 


12.8 


2.5-5.0 (3.8) 


11.2 


ii-n'"^ (14) 


1.4 


64.1 


5.8 


EXO0331+530/BQCam 


15.4 


7 


34.25 


20- 


1.4 


81.8 


5.4 


KS 1947+300 


14.2 


10.0 


40.4 


17.5'" 


1.4 


61.3 


4.5 


4U0115+634/V635 Cas 


16.3 


7-8 (7.5) 


24.3 


18' 


1.4 


58.8 


4.2 


IE 1145.1-6141/V830 Cen 


13.1 


8 


14.4 


10" 


1.4 


32.6 


4.0 


SAX J2103. 5+4545 


14.2 


6.5 


12.68 


20.0^° 


1.4 


45.4 


3.0 


4U 1907+09 


16.4 


5 


8.38 


26-27.9'^ (27) 


1.4 


49.2 


2.4 


4U 1538-52/QV Nor 


14.4 


4.5-6.4 (5.5) 


3.73 


19gl2.13 


1.4 


23.7 


1.6 


XTE J 1855-026 


~15 


10 


6.07 


2514,6 


1.4 


19.4 


1.0 


Cen X-3A'779 Cen 


13.3 


9 


2.09 


20.2* 


1.34' 


9.9 


0.6 



"The source distance given in lLiu, van Paradiis & van den Heuvdl 420061) and references therein. If[Li u, van Paradi Ts~& van den Heuvell 420061) give a range of 
values, then the value we adopt in this work is given in parentheses. If lLiu, va n Paradi is & v an den Heuvel (2006) give multiple values, then the value we adopt 
comes from the literature. Superscripts indicate the following references: (aa 1SteeIeTtan TT998); (hh] Kaper, van der Meer & Najarro 1 2006) 

^From iLiu, van Paradiis & van den Heuveli 12006.) and references therein. 

'^The estimate for the mass of the companion star from the literature. If a paper indicates a ran^ e of values, then the value we adopt in this work 
is giv en in parentheses. Superscripts indicate the following references: (1 Negueruela & Okazaki ('200 1"); (2'Negueruela et al.' f 1999); O^'Clark et alj 
<200Tb: (4lW ang & Gies (1998); (5; Kaaret et al. (1999); (6,Cox (2000); (7; Barziv et al. (2001); (8;van der Meer et al. (2007); (9; Ferrigno et al. (200^ 
(lO'Kaper, van der Meer & Najarro (2006); ( 11 Okazak i & Negueruela 12001.) ; ( 12 .Revnolds. Bell & Hilditch a992); (13.'Clark. 12000) ; (1 4jCorbet & MuiS 
002); (15 Cox, Kaper &Mokiem (2005); (16 VeiT ecchia et alj <2003) ~('17 1Wilson et alj 420031) ; (IS INegueruela et alj 420031) : d9 1CoeetalJ 419881) 
(20:iReig et al. ( 2004) 

''We assume a neutron s tar mass of I.4M(;) unless there is a value g ive n in the literature. Superscri pts indicate the following references: (w lBarziv et alj 4200 Ih ; 
(Ai 'iOuaintrell elal] 420031) ; (v IKaper. van der Meer & NaiaiTd 42006) ; (z lvan der Meer et alj UWH\ 



optical light comes from one of the binary components. The 
targets should be known to harbor NSs, and evidence for this 
can either come from the detection of X-ray pulsations or the 
presence of type I X-ray bursts. Finally, it is useful to know 
the orbital periods for planning the observations, and orbital 
periods as well as other parameters, such as the source dis- 
tance and the companion mass, must be constrained in order 
to estimate the angular sizes of the orbits. 

While one could consider High-Mass X-ray Binaries 
(HMXBs) and Low-Mass X-ray Binaries (LMXBs), it is clear 
that the former will provide the best opportunities. A large 
number of HMXBs are bright {V ^ 5-20), most of the optical 
light comes from the companion star, and their long orbital 
periods (days to a year) mean that the separation between the 
NS and the companion is large. In addition, many HMXBs 
are X-ray pulsars and for a large fraction of these the projected 
size of the NS orbit (a ^ sin /, where / is the binary inclination) 
has been measured, which, when combined with astrometric 
measurements, allows for a direct NS mass measurement. 

Thus, in selecting sources, we have started 
with the catalog of Galactic H MXB s compiled by 
iLiu. van Paradijs & van den Heuvell (l2006h . This cata- 
log includes 114 HMXBs, including 38 sources for which 
both orbital periods and NS spin periods have been measured. 
These HMXBs range from being very bright in the optical 
(V = 6) to not having known optical counterparts, and we 
further cut the list down to 27 sources by requiring that the 
sources have a V-magnitude brighter than 20. While all 
of these sources are known X-ray pulsars, o^ sin/ has not 
been measured in all cases. As this parameter is required 
for a direct NS mass measurement, we finally select the 
17 HMXBs that meet the above criteria and also have a 



measurement of the projected size of the NS orbit. 

The 17 sources are listed in Table [T] along with their V- 
band magnitudes and the information relevant to estimating 
the angular size of the companion's orbit. If the semi-major 
axis of the binary is a, then the corresponding angle subtended 
at the source distance d is p = tan"' (a/d), which is given by 

p = 2.35Mas4;,P'/'hr(^comp+MNs)'/' , (1) 

where t/^pc is the distance to the source in kpc, Porb.hr is 
the orbital period in hours, and M^omp and Mns are the 
masses of the companion and the NS, respectively, in units 
of Mq. Then, the companion's orbital signature is given by 
/Ocomp = pMns/(Mns +M cnmp)- In Table ffl the values of d 
and forb are taken from Liu, van Paradiis & van den Heuvell 
(2006). For systems where estimates of Mns and Mcomp have 
been given in the literature, that value is given along with 
the reference. Otherwise, companion masses are estimated 
using the spectral type (or range of spectr al types) given in 
iLiu. van Paradijs & van den Heuvell (l2006|) a nd the tables of 
stellar masses given in Chapter 15 oflCoxl (2000), and NS 
masses are assumed to be 1.4Mq. If some of these systems 
actually contain over-massive neutron stars, then pcomp will be 
larger than our estimates in these cases. 

3. SIMULATIONS 

3.L The Path of the HMXB Photocenter 

To represent motions of the target photocenters on the sky, 
we define a linear coordinate system such that x = Aq;cos(5 
and y = AS where a and & are, respectively, the Right As- 
cension and Declination of the sources at any time, t. The 
following equations represent the motion of the target photo- 
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centers 

= xo + (f - fo)M.i- + 7r cos sin (a© -ao)+ -^orbital and (2) 
y = yo + (f-fo)My+7rcos(5osinaocos(ao-Q;o)+}'oibitai , (3) 

where quantities with subscript "0" are evaluated at time to, 
Uq and ^0 are the R.A. and Decl. for the sun, and Xoibitai 
and joibitai represent the terms in the x- and y-directions due to 
orbital motion. Following Lindegren ( 1997) and section 2.3. 4 
of Volume 1 of the Hipparcos Catalog (.Ferryman et al.lfT997h . 
the orbital terms are 

x,,ut.i = BX(t) + GY(t) and y,,ut,i=AX(t) + FY(t) , (4) 

where A, B, F, and G are the Thiele-Innes elements (Gree^ 
[19851). which are a linearized alternative to the standard 
Campbell orbital elements, and depend on the angle sub- 
tended by the semi-major axis of the photocenter, pphot, the 
argument of periastron, w, the position angle of the node, 57, 
and the binary inclination, /. 
The time-dependent parameters are given by 

X(f) = cos£'-e and y(f) = (1 -e^)'''^sin£' , (5) 

where e is the orbital eccentricity, and E is the eccentric 
anomaly, which is given by 

27r 

— (f-r) = £'-esin£' , (6) 

where P is the orbital parameter, and T is the time of perias- 
tron passage. 

Thus, the path of the photocenter is determined by 12 pa- 
rameters: 2 position, 2 proper motion, 1 parallax, and 7 or- 
bital. As mentioned above, almost all of the light from an 
HMXB comes from the high-mass stellar companion, plac- 
ing the photocenter at the location of the companion star 
dCoughlin et al.ll20Tol) . In the following, we make the approx- 
imation that 100% of the light comes from the companion so 
that pphot corresponds to the semi-major axis of the compan- 
ion star itself (pcomp)- Figure [T] shows the path of the com- 
panion star for the HMXB X Per over a period of one year, 
illustrating the different contributions to the motion. 

3.2. SIM Lite Astrometry Measurements 

Here, we consider the specific capabilities of SIM Lite in 
order to carry out our simulations. SIM Lite will measure 
positions of grid stars across the sky, which will provide 
global astrometric measurements to a limiting accuracy of 
4 /ias and it will use reference stars around selected 

targets (such as the HMXBs) to provide a co-ordinate system 
covering a smaller region (~2°) where relative astrometric 
measurements will have a limiting accuracy of 1 /xas (l-rr), 
and these are referred to as "Narrow Angle" measurements 
dNemati & Mora les 2009: lUnwinl l2009). It is the relative un- 
certainty that is important for considering how accurately we 
will be able to measure the HMXB orbital parameters, and we 
use information from the SIM Lite project about the accuracy 
of these Narrow Angle measurements. However, parallax and 
proper motion measurements require global astrometry. 

In addition, for our simulations, we consider the fact that 
the SIM Lite measurements provide 1-dimensional positions. 
While the path of the companion star is fully defined by x and 
y, one must consider the baseline orientation for each mea- 
surement. We use the angle 6 to define the orientation of the 
baseline, where 6 is measured from the x-axis in the counter- 
clockwise direction. The quantity that is actually measured is 
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Fig. 1 . — The best estimate for the path of the companion star and pho- 
tocenter of the High-Mass X-ray Binary X Per over a period of 1 year. The 
dotted line shows the motion due only to the proper motion, and the dashed 
line shows the motion due to proper motion and parallax, assuming a distance 
of 1000 pc. The solid line includes these two effects as well as orbital motion. 
For X Per, the orbital period is 250 days, and we estimate that the semi-major 
axis of its orbit will subtend an angle of 165 ^tas. 



angle between the position of the target after it is projected 
onto the baseline and the origin of the x-y co-ordinate system, 
which is given by r = xcosO+ysmO. For the observing cam- 
paigns described below with many measurements, a baseline 
angle is chosen at random for each observation. 

Finally, for the simulations, we determine the uncertainty 
in each measurement using the on-line SIM Lite "Differential 
Astrometry Performance Estimator" (DAPE)"*. This tool al- 
lows the user to put in information about their target as well 
as details about the SIM Lite observations, and it determines 
a "single measurement accuracy" (SMA), which we use for 
our simulations, as well as the amount of time required for 
the observation. An example of the input and output from this 
tool is shown in Table|2] Once the SMA has been determined 
for a given target, a random number generator is used to se- 
lect a number from a Gaussian distribution with a equal to the 
SMA, and adding this number to r results in the measured an- 
gle r„,eas- Thus, each visit to the target generates the following 
information: r,„eas, the uncertainty in r,„eas, which is equal to 
the SMA, a timestamp (f), and a baseline angle iff). 

3.3. Fitting for Orbital Motion 

Once the simulated data are produced, they are fitted with 
a function that accounts for proper motions and parallaxes for 
the target and reference stars as well as orbital parameters for 
the target. The functional form is / = xcos6'+3'sin0, where 
X and y come from Equations 2 and 3 above. In total, the 
equation includes 1 1 free parameters: five are non-orbital (xq, 
ya, M.v, My, tt), and six are related to the binary orbit (pcomp, ^, 
Vl, i, e, and T). For the targets of interest, the orbital period is 
known very accurately, and we fix this parameter to the known 
value rather than leaving it as a free parameter. 

We wrote code in Interactive Data Language (IDL) to per- 
form the fits, and used the ID L routine curvefit. This routine 
uses the Marquardt method dBevington & RobinsonI [T992"). 
which incorporates a gradient-expansion algorithm and com- 
putes a least-squares fit. It is necessary to carry out the 
computations using double-precision arithmatic with strin- 

"* See |http://mscws4.ipac.caltech.edu/simtools_v2/| 
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gent convergence criteria. Although curvefit includes an op- 
tion to estimate rather than exactly compute the function 
derivatives for each iteration, we obtained better results when 
computing the exact derivatives. 

3.4. Cross-Check with Other Simulation Software 

We checked the software described above by comparing it 
to the code used for simulations of SIM Lite observations of 
planets (Traub et al. 2010). The planet code was developed as 
part of the SIM Double Blind Test and uses realistic cadences 
and baselines but assumes 2-dimensional observations. There 
were 5 teams involved in the Double Blind Test, including 
one team led by one of us (MWM). Table |3] shows that there 
is, in general, excellent agreement between the two sets of 
code (the planet code and the new code written for this work). 
Although there are a couple differences, they are easy to un- 
derstand. One difference is that the proper motions are sig- 
nificantly more accurate for the planet code, and this is due 
to the fact that the simulated observations are spread over 5 
years rather than 1 year. Second, the other parameters all have 
slightly smaller uncertainties with the planet code because the 
SMAs were smaller in a previous version of the DAPE calcu- 
lator 

4. RESULTS 

4. 1 . HMXBs for Simulations and Observing Parameters 

Figure |2] shows the predicted orbital astrometric sig- 
natures and the V-band magnitudes for the 17 HMXBs 
listed in Table [T] comparing their location in the plot to 
the approximate sen sitivity limit for SIM Lite taken from 
iTomsick. Shaklan & Pani 62009 ). The limit is based on 
40 hours of mission time and is conservative in the sense that 
it corresponds to a signal-to-noise ratio of 10. Thus, while 
some of the sources that are close to the limit may still have 
orbits that are measurable by SIM Lite, in this work, we fo- 
cus on the six systems that are clearly above the line: X Per, 
V725 Tau, Vela X-1, GX 301-2, SAX J0635.2H-0533, and 
V830 Cen. 

For each of our targets, we used the DAPE tool (see above) 
to determine a possible plan for observing campaigns that 
would use 48 hours of mission time (per target). We used 



Table 3 

Comparison of Measurement Uncertainties for 
Vela X-1 



Planet Code New Code 

Pai'ameter Description Units Uncertainty" Uncertainty* 



e 


Eccentricity 




±0.04 


±0.041 


i 


Inclination 


degrees 


±1.14 


±1.29 


w 


Argument of Periastron 


degrees 


±20 


±27 


to 


Time of Periastron Passage 


days 


±0.5 


±0.67 


n 


Position angle of the node 


degrees 


±1 


±1.4 


pcomp 


Semi-major axis 


^as 


±0.2 


±0.22 


d 


Distance 


pc 


±0.5'- 


±0.7^ 




Reference Position 


arcsec 


±4 X lO-'' 


±5.4 X 10- 


yo 


Reference Position 


arcsec 


±4 X 10-' 


±5.4 X 10" 


V.v 


Proper Motion 


arcsec/yr ±1.1 X 10 ' 


±5.2 X 10- 


Vv 


Proper Motion 


arcsec/yr ±1.1 X 10"' 


±5.8 X 10- 



"These are the 68% confidence (I-ct) uncertainties on the Vela X-1 parameters calcu- 
lated with the planet code that was written for the SIM Double Blind Test. These simula- 
tions include 100 observations made over 5 years. The observations are 2-dimensional, 
and the single measurement accuracy in each dimension is 1.27 /^as. 

^These are the 68% confidence uncertainties on the Vela X-1 parameters calculated 
with the new code written for this work. These simulations include 200 observations 
made over 1 year. The observations are 1 -dimensional, and the single measurement 
accuracy is 1.47 /xas. 

^Although this is the error from the fit to the NaiTow Angle observations, the actual 
uncertainty on the distance will depend on the tie-in to the absolute reference frame. 



the parameters shown in Table |2] except for the science target 
V-magnitude, which depends on the target, and the science 
target integration time per chop''. For the brightest four of 
our six targets, we were able to obtain low SMAs (<1.8 /ias) 
with the integration time set to 30 s (see Table lU. For the 
fainter two targets, we increased the integration time to 60 s, 
giving SMAs of 2.39 and 2.66 ^^as for SAX J0635.2H-0533 
and V830 Cen, respectively. The trade-off is having smaller 
numbers of observations; however, because we still use 30 s 
for the reference target integration time per chop and about a 
third of the time for the campaign is overhead (see Table |2]i, 
we only need to decrease the number of visits in the observing 
campaign from 120 to 90 for the fainter two targets. 
While detection of orbital motion may still be possible for 



Table 2 

Differential Astrometry Performance Estimator for 
SIM Lite: INPUT AND OUTPUT 



Quantity 


Value 


Input parameters 


Science target V-magnitude 


7 


Science target B-V color index 


0.4 


Characteristic reference target V-magnitude 


10 


Characteristic reference target B-V color index 


0.4 


Target-reference sky separation 


1° 


Number of reference targets 


4 


Number of chop cycles 


4 


Science target integration time per chop 


30 s 


Reference target integration time per chop 


30 s 


Output parameters 


Single measurement accuracy 


1.47 /ias 


Science target integration time per visit 


480 s 


Reference target integration time per visit 


480 s 


Overhead time per visit 


480 s 


Total mission time per visit 


1,440 s 



Table 4 

Observing Parameters for a Campaign with 48 Hours 
OF Mission Time 







Science Target 


Single 








Integration Time 


Measurement 


Number of 


Source Name 


V 


Per Chop (s)" 


Accuracy (/ias)'^ 


Observations'^ 


4U 0352-^30/X Per 


6.6 


30 


1.45 


120 


3A 0535-^262/V725 Tau 


9.6 


30 


1.56 


120 


Vela X-l/GP Vel 


6.9 


30 


1.47 


120 


GX 301-2/BPCr-u 


10.8 


30 


1.76 


120 


SAXJ0635.2-^0533 


12.8 


60 


2.39 


90 


IE 1145.1-6141/V830 Cen 


13.1 


60 


2.66 


90 



"This value was varied according to the V-magnitude of the source, but the 
other input parameters were left at the values shown in Table|2] 

''These values come from the Differential Astrometry Performance Estima- 
tor for SIM Lite 

^In each case, the number of observations was adjusted to obtain a total 
mission time, including overheads, of 48 hours. 

^ In making a measurement of an angle between sources, the SIM Lite 
instrument "chops" back-and-forth betw een the target and the reference star 
by changing the length of the delay line jUnwinll2009l) . 
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Fig. 2. — Predicted astrometric signatures from orbital motion vs. V-band magnitude for 17 neutron star HMXBs for which the projected size of the neutron 
star's orbit (Oy sin i) has been measured via X-ray pulsations. The dashed line shows the threshold for detection of orbital motion in 40 tiours of SIM Lite mission 
time. As described in iTomsick, Shaklan & Pan 2009), the threshold is defined as the level at which the system's semi-major axis is 10 times larger than the 
astrometry noise per observation (i.e., the single-measurement accuracy) divided by the square root of the number of observations. 



several of the other HMXBs close to the sensitivity line in 
FigurelH such a detection is not feasible in the extreme cases. 
For a source as faint as EXO 2030+375 (V = 19.7), even if 
one increases the science target integration time per chop to 
2 hours (the maximum allowed by the DAPE tool), the SMA 
is 10 /ias, which is larger than the 7 fias size of the orbit (see 
Table[Tli. Although one could conceivably still detect the orbit 
with a large number of visits, even as few as 20 visits would 
require nearly 1 month of mission time. Thus, detecting the 
orbit would likely take several months and would clearly not 
be a good use of SIM Lite time. 

4.2. Fitting Results 

For each of the six HMXBs, we simulated 2,000 SIM Lite 
observing campaigns and fitted the data with the 1 1 parame- 
ter model described above (note that there are 1 1 rather than 
12 free parameters because we always fixed the orbital period 
to the known value). For each HMXB and each parameter, 
we produced histograms of the 2,000 best fit values obtained. 
For all of the parameters of all of the HMXBs except for 
V830 Cen, the histograms appear to be normally distributed. 



allowing us to fit each histogram with a Gaussian and use the 
a value to estimate the uncertainties on each parameter that 
SIM Lite will be able to obtain. For V830 Cen, the histogram 
for the eccentricity parameter is not normally distributed. 
We suspect that the parameter is not well-constrained be- 
cause V830 Cen has the smallest estimated semi-major axis, 
Pcomp = 4.0 /ias, and we discuss the case of V830 Cen further 
below. 

We report the simulation results for the other five HMXBs 
in Table |5] including the six orbital parameters as well as the 
parallax measurement. SIM Lite will greatly improve con- 
straints on the binary inclinations, /. While the current uncer- 
tainties are on the order of ~10°, with a large systematic com- 
ponent, SIM Lite will constrain / to within ±0.28°, ±0.98°, 
±1.58°, and ±2.50° for X Per, V725 Tau, Vela X-1, and 
GX 301-2, respectively. Also, SIM Lite will provide the first 
direct measurements of the angular size of the companion's 
semi-major axis, pcomp- For example, for X Per and Vela X-1, 
our simulations indicate measurements of 165 ±0.45 /ias and 
10.2 ±0.28 /ias. 

For each HMXB, Table |5] provides a measurement of the 
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parallax, tt, and then an actual measurement of the distance. 
The uncertainty shown with tt is that obtained from the fit to 
the data. However, there is also an uncertainty from the fact 
that the reference frame will only be known to within 4 /ias, 
and both error components are given with the distance mea- 
surement. If X Per is at a distance of 1000 pc, we predict that 
SIM Lite will measure its distance to an accuracy of ^5 pc. 

For Vela X-1 and GX 301 -2, the X-ray pulsation measure- 
ments (Bil dsten et al.l [19971) give values of e and lo that are 
much more accurate than our simulations predict for SIM Lite. 
Thus, we have re-fit our simulated data for these two systems 
and also for the two fainter sources, SAX J0635.2H-0533 and 
V830 Cen, keeping e and a; as fixed parameters. Once this 
is done, the V830 Cen histograms appear to be normally dis- 
tributed. Table |6] shows the fitting results for these four sys- 
tems. The improvements in the constraints on the most im- 
portant parameters (/, Pcomp, and d) are relatively small. The 
largest improvement is in the measurement of the time of pe- 
riastron passage (fo), which is expected since this parameter is 
related to e. 

4.3. Direct Neutron Star Mass Measurements 

It is possible to write an expression for the mass of a neutron 
star, Mns, in an HMXB in terms of five directly measurable 
quantities. Two of the quantities: forb and the projected linear 
size of the neutron star orbit (Of sin/); are measured by instru- 
ments other than S'/ML/fe. The other three quantities: pcomp, 
/, and d\ will be measured by SIM Lite as described above. 
Starting from the standard orbital equations that are typically 
used for measurements of the components of binary systems 
(e.g.. Equation 5.3 in .Charles & Co&.2006.) . we derive the fol- 
lowing equation 

47r^ J tanpcompr, • •^ , j , ■ -^l /^^ 

Mns = -p^ —i — - sm i) + d tan p^omp sm ;] ^ , (7) 

(j^orb sm ; 

where Mjvs, Porh, d, and axsini are in the same units as the 
gravitational constant, G, e.g., CGS units, and / and pcomp 
are angles. For X-ray pulsars, the projected size of the neu- 
tron star orbit can often be determined to very high accu- 
racy, e.g., axSini = 1 13.89 ±0.13 light-seconds for Vela X-1 
dBildsten et al.|[T997 l). By obtaining long (years) time base- 
line X-ray or optical observations. X-ray pulsar orbital peri- 
ods are, for our purposes, known with negligible uncertain- 
ties. 

In addition to the SIM L/fe-measured parameters in the top 
parts of Table |5] the values of a j sin / are given for the 5 
HMXB systems, and this, in turn, is used to calculate Mffs- 
For X Per, our simulations predict that SIM Lite will measure 
the neutron star mass to ~2.5% (M^s = 1.40 ±0.035 Mq, 
where the uncertainty corresponds to the 68% (1-cr) confi- 
dence level). The second-best measurement would be ob- 
tained for Vela X-1, for which we predict a ^6.5% uncertainty 
on Mns- It should be noted that these values consider both 
components (from the fit and from the reference frame) of the 
distance uncertainty. For each HMXB, the contribution from 
each parameter (d, Pcomp, and a^sin/) to the uncertainty on 
Mns is given in Tables \5\ and |6] For V725 Tau, the largest 
contribution to the uncertainty comes from the a^sini term, 
so that the future measurement of M^s can be substantially 
improved with X-ray observations of this system. 

5. DISCUSSION AND CONCLUSIONS 

The results of these simulations show that SIM Lite will 
provide excellent constraints on the masses of neutron stars 



in HMXBs. There is currently no direct method for obtaining 
the binary inclinations (/) of HMXBs, and this will be a major 
improvement in the measurements. This is especially true for 
non-eclipsing systems like X Per and V725 Tau, for which / 
is thought to be near 30°, but it is currently only estimated to 
about ±10-15°. With such a large uncertainty in /, we do not 
have any current estimate of the NS masses in these systems, 
which is why we assume 1.4M0 for the simulations. Since 
these are both accreting systems, there is a strong possibility 
that the NS masses are significantly higher, which could allow 
for a constraint on the NS EOS. 

The current constraints on / are better for Vela X- 1 , which 
is an eclipsing system, but HMXBs with relatively short or- 
bital periods can be eclipsing for a wide range of inclinations; 
thus, the improvement in the measurement of / that SIM Lite 
will provide is important. There are akeady suggestions that 
the Vela X- 1 NS is over-massive, and if SIM Lite finds that the 
NS mass is, e.g., 2.00±0.13 M0, as fo und with our simula- 
tions, this would rule out many NS EOSs (iLattimer & PrakasE 
I2OOI . 

There are also possibilities for further improvements to the 
accuracy of the NS mass measurements estimated in Tables |5] 
and |6] We assume in our simulations that we will use SIM 
Lite to determine the times of periastron (fo), but if contempo- 
raneous X-ray observations can be made, it will be possible 
to accurately determine this time using the X-ray pulsations. 
Also, we find that, after SIM Lite measurements of V725 Tau, 
the uncertainty in M^s will be dominated by the measurement 
of fl V sin /; thus. X-ray measurements to improve the measure- 
ment of this parameter could ultimately lead to a mass con- 
straint that is nearly as good as that for X Per Finally, one 
more piece of information that we have not considered in this 
work is the radial velocity semi-amplitude of the companion 
star (Kcomp)- This parameter is related to the projected size 
of the orbit, and a cross-check between this parameter and 
the astrometric parameters (pcomp, d, i, and e) will be possible 
(e.g., Pour baix & Jorissen 2000). 

Although we focus on obtaining NS masses via orbital mea- 
surements in this work, it should also be pointed out that SIM 
Lite can also contribute to constraining EOSs by improving 
measurements of NS radii. Two techniques for measuring the 
radii of NSs in Low-Mass X-ray Binaries use measurements 
of the X-ray emission from the neutron star surface: one is to 
observe the thermal X-rays from the surface of a NS whe n 
the accretion rate is very low (e.g., iRutledge e t al.l l2002h ; 
and another is to measure the thermal emissi on near the 
end of a type I X-ray burs t (iGalloway. Ozel & PsaltisI 120081: 
lOzel. Baym & Guveil2010l) . In both cases, the distance to the 
LMXB is the largest uncertainty in the NS radius measure- 
ment (see the references above and iTomsick. Shaklan & PanI 
120091) . and SIM Lite will provide accurate distances for a large 
number of LMXBs. 

Finally, it is important to note that our requirement in 
this work that a^sini be a measured parameter eliminates 
many interest ing HMXBs that will be excel lent targets for 
SIM Lite. In iTomsick. Shaklan & PanI (|2009|) . we find ~20 
HMXBs for which the predicted astrometric signatures will 
be large enough for orbital measurements. This includes 
we ll-known sources such as the black hole system Cyg X- 
1 (ICaballero-Nie ves et al.' '200 9|) and the Ukely black h ole 
system SS 433 (^lundell. Bowl er & Schmidtobreickl2008l) as 
well as systems where it is not clear whether the compact ob- 
ject is a black hole or a NS. These include the interesting case 
of 4U 1700-377, which is thought to be a NS based on its X- 
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Table 5 

Simulation Results and Mass Measurements (all parameters free) 



Source Name 


XPer 


V725 Tau 


Vela X-1 


GX 301-2 


SAX J0635. 2+0533 


SIM Lite orbital measurements 


e 


0.1 1 ±0.0026 


0.47 ± 0.0099 


0.0898 ± 0.050 


0.462 ±0.068 


0.29±0.15 




29.82 ±0.28 


30.50 ±0.98 


82.52± 1.58 


67.40 ±2.50 


49. 12 ±7.84 




108 ± 1.5 


310 ±2.2 


332.59 ±33.6 


130.4±8.2 


356 ±36 


to (days) 


±0.99 


±0.33 


±0.84 


±0.84 


±0.88 


no 


45 ±0.5 


45 ± 2.0 


45 ± 1.6 


45 ± 2.7 


45± 12 


Pcomp (M^IS) 


165 ±0.45 


49.0 ±0.41 


10.2 ±0.28 


9.8 ±0.49 


5.8 ±0.56 


SIM Lite distance measurement 


TT (//as) 


1000 ±0.87 


500 ± 0.37 


526 ±0.28 


278 ±0.31 


263 ±0.68 


rf(pc) 


1000(±0.87)(±4) 


2000(±1.5)(±16) 


I900(±1.0)(±14) 


3500(±4.0)(±52) 


3800(±9.9)(±58) 


X-ray measurement of the projected size of the pulsar's orbit 


fly sin / (It-s)'' 


454.0 ±4.0' 


267 ± 13' 


113.89±0.I3- 


368.3 ±3.7^ 


83± ll' 



Derived neutron star mass measurement and enor contributions 



Mns (Mq) 


1.40 ±0.035 


1.40±0.15 


2.00 ±0.1 3 


1.85 ±0.21 


1.40 ±0.59 


SMNS,di„a:,ce (Mg) 


0.008 


0.014 


0.018 


0.032 


0.029 


l5M«S,pcomp (M©) 


0.004 


0.014 


0.062 


0.103 


0.156 


SMnsj (Mq) 


0.025 


0.087 


0.110 


0.176 


0.463 




0.023 


0.125 


0.004 


0.036 


0.329 



"These are from the following: (lIDelgado-Mani etalj <200lh: (21Bildsten et alj (199%: OIKaaret. Cusumano & Saccol <2000l) 



Table 6 

Simulation Results and Mass Measurements (e and u} fixed) 



Source Name 


Vela X-1 


GX 301-2 


SAXJ0635.2-^0533 


V830 Cen 


SIM Lite orbital measurements 


e 


0.0898° 


0.462° 


0.29° 


0.20° 


iC) 


82.52 ± 1.55 


67.40 ±2.29 


49.12 ±7.78 


60.68 ± 9.87 




332.59° 


130.4° 


356° 


128° 


to (days) 


±0.040 


±0.27 


±0.31 


±0.47 


no 


45 ± 1.6 


45 ± 2.4 


45± 10 


45± 13 


Pcomp (^as) 


10.2 ±0.27 


9.8 ±0.37 


5.8 ±0.52 


4.0 ±0.56 


SIM Lite distance measurement 


TT (pas) 


526 ±0.27 


278 ±0.31 


263 ± 0.67 


125 ±0.51 


d(pc) 


1900(±1.0)(±14) 3500(±4.0)(±52) 


3800(±9.7)(±58) 


8000(±33)(±256) 


Derived neutron star mass measurement and enor contributions 






2.00 ±0.1 3 


1.85 ±0.19 


1.40 ±0.59 


1.40 ±0.57 


SMNS,dis,ancr (Mq 


) 0.018 


0.032 


0.029 


0.062 


*'>^«S,p„mp (Mq) 


1 0.061 


0.078 


0.145 


0.242 


SMNS,i (Mq) 


0.108 


0.161 


0.463 


0.516 


SMfis.nsini (Mq) 


0.004 


0.036 


0.329 


0.044 



"Fixed. 

*The value of Uxsini used for V830 Cen is 99.4 ±1.8 It-s jRav & Chakrabart^l2002l) . The values used for the other sources are the same as in Table|5] 



ray properties, but which has a compact object mass estimated 
at 2.4 M0. Also, precise compact object masses have not been 
determin ed for the gamma-ray binaries LS 1-1-61° 303 and 
LS 5039 JPubu s 2006), but given that there are only a small 
number of high-mass binaries known to produce gamma-ray 
emission, determining the type of compact object in these sys- 
tems is of great interest. 
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